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Abstract. We present in this paper a catalog of reference stars suitable for calibrating infrared interferometric observations. 
In the K band, visibilities can be calibrated with a precision of 1 % on baselines up to 200 meters for the whole sky, and up 
to 300 meters for some part of the sky. This work, extending to longer baselines a previous catalog compiled by Borde et 
al. (2002), is particularly well adapted to hectometric-class interferometers such as the Very Large Telescope Interferometer 
(VLTI, Glindemann et al. 2003) or the CHARA array (ten Brummelaar et al. 2003) when observing well resolved, high surface 
brightness objects (A^ S 8). We use the absolute spectro-photometric calibration method introduced by Cohen et al. (1999) to 
derive the angular diameters of our new set of 948 G8-M0 calibrator stars extracted from IRAS, 2MASS and MSX catalogs. 
Angular stellar diameters range from 0.6 mas to l.Smas (median is 1.1 mas) with a median precision of 1.35 %. For both the 
northern and southern hemispheres, the closest calibrator star is always less than 10° away. 
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1. Introduction 

Long baseline stellar interferometers measure the amount of 
coherence (ie the squared modulus of the coherence factor, 
sometimes also called the raw visibility amplitude, or raw 
fringe contrast) between any pair of their subapertures. This 
quantity needs to be calibrated in order to yield the true visibil- 
ity of the source, which is the modulus of the Fourier transform 
of the object's intensity distribution at the spatial frequency 
B/A defined by the projected baseline B and the wavelength 
A. The squared visibility V~ of an object is derived from the 
measure fr of its coherence factor, and from an interferometric 
efficiency factor (also called transfer function), which ac- 
counts for the coherence losses caused by imperfections of the 
instrument and by the Earth's turbulent atmosphere; 



owing to changing instrumental and atmospheric conditions, 
it has to be frequently sampled, and then interpolated for the 
time of the science target observation (one can see for example 
Perrin (1998) for details of an interpolation strategy) to ensure 
a reliable calibration.The determination is done by observ- 
ing a reference source (yu^^j ) with a known squared visibility 
(y2^, ), for which 



rj-2 ^ ^cal. 



(2) 



r2 



(1) 



From Eq. 1, it appears that the errors on and fi^ con- 
tribute equally to the accuracy of the object squared visibil- 
ity V~. As the interferometric efficiency varies during the night 
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In practice, observations of the science target are inter- 
leaved with observations of reference stars, referred to as cali- 
brator stars or calibrators for short. The choice of a calibrator 
is critical as one should be able to predict its intrinsic visibility 
at the projected baseline of the interferometer, with an accu- 
racy which should be at least as good as the measurement of 
the raw visibility (more about this in Section 2). For this pur- 
pose, Borde et al. (2002) carefully extracted a catalog of cali- 
brator stars from the spectro-photometric reference stellar net- 
work compiled by Cohen et al. (1999), hereafter C99. Because 
of the estimated angular diameter of the selected stars (typically 
2.3 mas with an error of 1.2 %), this catalog, hereafter Cat. 1, 
provides calibrators for interferometric baselines up to ~ 100 m 
in the near infrared. The advent of large interferometers, such 
as the VLTI and CHARA arrays with maximum basehnes of 
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202 and 330 m, respectively, motivated us to compile a catalog 
for longer baselines, hereafter Cat. 2. 

In Sect. 2, we explain why the sensitivity of current inter- 
ferometers leads us to take as calibrators partly resolved stars 
instead of point-like sources, in the bright star regime. Then, we 
argue that these stars can be safely modeled as uniform disks 
provided they are carefully selected, and we review two differ- 
ent strategies to make this selection. Last, we explain how to 
decide if a calibrator is appropriate for a given observation. In 
Sect. 4, we describe how we built Cat. 2 by using C99's method 
and infrared photometric measurements for stars fainter than in 
Cat. 1. Finally, in Sect. 5, we discuss the main characteristics of 
Cat. 2. In particular we demonstrate the relevance of this new 
catalog for the VLTI/VISA (the VLT interferometer sub-array, 
involving the 1.8 m auxiliary telescopes) and the CHARA ar- 
ray. 

2. Choosing and modeling calibrators 

The ideal calibrator is a source for which the intrinsic visibility, 
using the same observational setup as the science target, can be 
perfectly predicted. Thus every calibrator needs to be described 
by a morphological model, so that its true squared visibility can 
be predicted for the full range of spatial frequencies that can be 
addressed by the interferometer. As reference sources are cho- 
sen among stellar objects, the most common models employed 
to describe a calibrator are, in increasing number of free pa- 
rameters: 

1. The point source (no free parameters): in that case V^^j = 1 
at any baseline; 

2. The uniform disk (UD model), with the angular UD diame- 
ter as the only free parameter: the monochromatic squared 
visibility is then given by the relationship: 



where /„ denotes the nth Bessel function of the first kind 
and X the spatial frequency, namely nB/A (B, baseline and 
A the observing wavelength); 
3. The limb darkened disk (LD model), with two or more free 
parameters. 

For simplicity and reliability, more complex calibra- 
tor models should be avoided. These include non centro- 
symmetric or variable morphologies like pulsating or flare 
stars, fast rotators, and binary stars. Indeed, a good calibrator 
needs to be: 

1 . Well modeled, in the sense that the model which is used to 
describe the morphology of the object is as simple as pos- 
sible, yet appropriate within the level of accuracy required 
for the visibility; 

2. Well known: the free parameter(s) in the model need to be 
known with sufficient accuracy; 

3. Stable: no time dependency of the morphology; 

4. Observable with the same setup and in the same conditions 
as the science target, e.g. close to the target on the sky, with 
comparable magnitude and preferably spectral type. 



It should be noted that these requirements are conflicting 
since in general they call for references much smaller than the 
science targets, yet having approximately the same color and 
apparent brightness, which is impossible for thermal sources. 
Therefore the choice of a calibrator will always be the result of 
a compromise. The first and last requirements exclude the point 
source model in almost every circumstance for long baselines, 
as it would imply an unrealistically small (hence faint) refer- 
ence. For example, in order to induce a bias smaller than 1 % 
{V^ > 99 %), a reference star considered as a point source with 
a 300 m baseline alA- 2.2 jum needs to have an angular diam- 
eter smaller than 0.1 mas. If the source is a KO giant, this corre- 
sponds to an infrared magnitude ^ 8, beyond the sensitivity 
limit of interferometers with small or medium size telescopes 
(less than 2 m in diameter) without fringe tracking devices. For 
interferometers using larger telescopes, with a sensitivity limit 
close to or better than K - 10, it appears possible to find real 
non-resolved sources (V^ > 99%), because sensitivity is no 
longer an issue. Still, one has to be careful for the following 
reasons. A A' = 10 star, earlier than M5, will definitely be un- 
resolved for a 200-m baseline interferometer operating in the 
near infrared, but one has to ensure that this particular candi- 
date follows the chosen model, i.e. a point source. Following 
our previous discussion, the well known item is out of concern 
for an unresolved star, in the sense that no parameter drives 
the model, but the well modeled one still applies. For exam- 
ple, the candidate may be binary, just like the majority of main 
sequence stars (eg see Quist and Lindegren 2000). Therefore, 
non-sensitivity limited interferometers still need well mod- 
eled (thus trusty) calibrators. Considering that faint candidates 
{K ^ 8) will arise from infrared photometric surveys, such as 
DENIS or 2MASS, for which very few supplementary data are 
available, the potential binarity will be very difficult to detect. 
It appears that there are two different cases for modern hecto- 
metric baseline stellar interferometers: the bright star regime 
{K $ 8) and the faint star regime {K ^ 8). Because one wants 
the interferometer to operate in the same regime for both the 
scientific target and the calibrator (eg for signal to noise ratio 
issue), we advocate that an extension to the Borde et al. (2002) 
catalog, containing only well modeled (uniform disk ensured) 
and well known (high accuracy estimation of the angular di- 
ameter), will not only be useful for medium aperture arrays op- 
erating in the near infrared, such as the CHARA Array or the 
Palomar Testbed Interferometer (Colavita et al. 1999), but also 
for the VLTI/VISA interferometer using AMBER in the bright 
star regime. 

It happens that quiet single stars can easily and correctly be 
modeled as uniform disks when the squared visibility exceeds 
Si 40 %, as the difference between a UD and a LD model is then 
smaller than 0. 1 % (see Fig. 1). Steflar angular diameters can be 
known from direct high-angular resolution measurements (in- 
terferometry or Lunar occultations) or estimated indirectly by 
(spectro-)photometry. We argue that a spectrophotometric esti- 
mate of the diameter is more suitable for reference purposes, 
because direct diameter measurements have only been per- 
formed on a limited number of sources (listed in the Catalog of 
High Angular Resolution Measurements (Richichi et al. 2004), 
hereafter CHARM), and most of these sources do not have the 
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Fig. 1. Relative difference between UD and LD squared visibilities 
{6V^/V~) with respect to spatial frequency times angular diameter. 
The UD diameter is adjusted to minimize the quadratic deviation be- 
tween UD and LD squared visibilities in the > 40 % domain. This 
particular LD model represents the most limb-darkened disks of our 
catalog, namely for MOIII stars in the J band. In the domain corre- 
sponding to > 40 %, the relative difference remains under 0. 1 %. It 
is even smaller for less limb-darkened stars. 



properties required for being a good calibrator at long base- 
lines. Besides, this interferometric data set is in essence hetero- 
geneous and more specifically, UD values are published usually 
at a single wavelength. Because of limb darkening, these val- 
ues are not readily valid for other band-passes. Therefore, in 
absence of an exhaustive and homogeneous set of direct diam- 
eter measurements of suitable calibrator stars (a certainly much 
needed, however intensive, observational program), we prefer 
in the remainder of this paper to apply an indirect method to 
determine the diameters for a set of stars that are selected so 
that they can be trusted as single and stable. 

At this point, two strategies have been proposed: 

1. A dynamic or open list approach: a software that queries 
existing databases for calibrator candidates in a given field 
is designed and distributed. The software is expected to 
help the observer in the down-selection process by using 
the information available in the databases. Then, it should 
provide either a direct estimate of the visibility at the time 
of the observation, or the parameter value(s) for a given 
visibility model (eg the UD angular diameter). Two exam- 
ples of this approach are Search Calibrator integrated in 
the software ASPRO' issued by the Jean-Marie Mariotti 
Center (Bonneau et al. (2004)), and getCaf issued by the 
Michelson Science Center; 



' http ://mariotti .uj f-grenoble. f r/~ aspro/ 
^ http://msc.caltech.edu/software/getCal/ 



2. A static or closed list approach: a catalog of carefully 
selected potential calibrators and including all the neces- 
sary information to compute the visibility for any rele- 
vant baseline and bandpass, is compiled and distributed. 
Such a list, the ESO Calvin TooP, has been made avail- 
able for the VLTI Mid-Infrared interferometer (MIDI, see 
Leinert et al. 2004). 

We prefer the closed list approach for three reasons: 

1. It enables a dedicated check on all of its entries, which 
would not be possible with a generic selection process; 

2. The catalog can be revised and updated as some of its en- 
tries are found - either serendipitously or by means of a 
systematic interferometric survey - to be unreliable as a 
reference (the most likely reason being that they are hith- 
erto unknown binaries), with the objective that eventually 
all bad entries will be removed; 

3. Selecting a calibrator is easier for the beginner in interfer- 
ometry since it can be hand-picked from a source list. 

However, the dynamic approach has its merits as it could 
provide either specific calibrators missing in our catalog (eg 
stars earlier than KO), or standard calibrators for the faint star 
regime (K ^ 8) where a closed list is unmanageable. 

3. The calibrator decision diagram 

The purpose of the calibrator decision diagram described here 
is to find the right calibrator for a given observation. In order to 
achieve the desired precision for a given instrumental config- 
uration, calibrators have also to be selected according to their 
diameters and diameter errors. For the UD model, the stellar 
angular diameter contains all the knowledge of the source, thus 
the squared visibility error associated with the calibrator error 



^ "^cal. > 



^cal. 



= 2x61- 



(4) 



As the knowledge of the calibrator should not be the lim- 
iting factor on the precision of the visibility measurement, we 
require for the calibrator choice that the relative error on the 
squared visibility, o-yi /y-^j , due to the uncertainty on the cali- 
brator diameter alone, <to/6, should not exceed the internal (in- 
strumental) error, cr^2 /yu^, that is to say 



I y2 



(5) 



limit 



For further discussions, we set the upper limit on inequal- 
ity (5) to 2 % as it implies on the relative error on the vis- 
ibility an upper limit of o-y/V - 1 %, the standard preci- 
sion benchmark for single-mode interferometers in the litera- 
ture. Examples of scientific results achieved with such a pre- 
cision include the study of the internal structure of Sirius by 
Kervella et al. (2003) or the test of the model atmosphere of a 
M4 giant by Wittkowski et al. (2004). 



http://www.eso.org/observing/etc/preview.html 



4 Merand et al.: A catalog of bri] 




0.1 1.0 10.0 



ae/e (%) 

Fig. 2. Calibrator decision diagram (logarithmic scales) fixing the 
wavelength to 2.2 yum and the baseline to 200 m: contours of the in- 
strumental precision as a function of the diameter and diameter error 
of the calibrator Stars below a given contour could be used at the la- 
beled precision. Vertical dashed lines delimit the area corresponding 
to the typical diameter error obtained through (spectro-)photometric 
determinations (1-5 %). 



Inequality 5 can be turned into a decision diagram that 
helps decide whether a calibrator is suitable or not for a 
given observation. In the framework of the UD model and 
leaving aside magnitude and spectral type considerations, 
a given source can be represented as a point in a diame- 
ter/diameter error plane, ie {0, erg Iff). We refer to this plane 
as the calibrator plane. As for the instrumental configuration, 
all needed information is embedded in the pair spatial fre- 
quency/instrumental precision, i.e. {x, (cr^2/yu^)iiniit) or alterna- 
tively in the triplet wavelength/projected baseline/instrumental 
precision, ie {A, B, {cr^i I ir^^ii). By holding two of these three 
quantities fixed, one can plot the third one as a function of the 
diameter and diameter error We call this plot the calibrator 
decision diagram. 

As a first example, we set on Fig. 2 the wavelength to 
2.2 yum, the baseline to 200 m, and we plot contours of the 
instrumental precision in the calibrator plane. In that case, if 
one wants a-yi /V^^j to be at most 2 %, the diagram shows that 
either < 1.3 mas with a 1 % error or < 0.65 mas with a 
5 % error is required. It is noteworthy that for a-yi jV^^y < 2 % 
and crg/9 > 1 %, the squared visibility stays above 40 % (corre- 
sponding to a 6* = 1.3 mas star, observed at /I = 2.2 jum, with a 
200 m baseline), which assures of the validity of the UD model 
as discussed previously. For this reason, we keep in the follow- 
ing these requirements on the squared visibility and angular 
diameter precisions. 
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Fig. 3. Calibrator decision diagram (logarithmic scales) fixing the 
wavelength to 2.2 ^m and the instrumental precision to 2 %: contours 
of the baseline as a function of the diameter and diameter error of the 
calibrator Stars below a given contour can be used at the labeled base- 
line. Stars from Cat. 1 appear as gray squares and stars from CHARM 
as black dots. The double arrow indicates the direction we chose to 
follow to expand Cat. 1. The 50, 100, 200 and 300 m curves in the K 
band (A = 2.2/j/n) corresponds respectively to 35, 75, 145 and 220m 
in the H band (A = 1.6/jm) and 25, 55, 110 and 165m in the J band 
(A = 1.2/im). 



As a second example, we set in Fig. 3 the wavelength to 
2.2 /im, the instrumental precision to 2 %, and we plot con- 
tours of the baseline in the calibrator plane. Black dots repre- 
sent stars flagged as single and with no shell in CHARM, and 
gray squares represent stars from Cat. 1 . Stars located below a 
given contour can be used as calibrators at the corresponding 
baseline. It appears from this plot that CHARM provides cal- 
ibrators for projected baselines up to ~ 70 m, and Cat. 1 for 
baselines up to ~ 100 m, not being too stringent on the sky 
coverage. 

4. Building the catalog 

4.1. How to expand Cat. 1: better models or smaller 
stars? 

Looking at Fig. 3, there are two ways in which Cat. 1 could be 
expanded to satisfy the 2 % precision at 200 m baselines: 

1. Decrease the diameter errors by improving the calibrator 
models (going leftwards on the diagram); 

2. Decrease the stellar diameters by finding smaller stars (go- 
ing downwards on the diagram). 

The first solution necessitates working with more resolved 
stars which necessitates more elaborate models than the UD, 
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along with high-quahty observations to determine the parame- 
ters of these models. Keeping with our objectives of simplicity 
and reliability, we favor the second approach, which consists 
in going only to smaller stars since we guaranty hereby the va- 
lidity of the UD model. On the other hand, we have to be very 
careful not to degrade the error. The second solution features 
the drawback of working with fainter stars. Indeed, the angular 
diameter logarithmic scale on Fig. 3, can be viewed as a linear 
magnitude scale for a given spectral type (a given temperature). 

4.2. Pushing the spectro-photometric method further 

Through their whole series of papers Cohen et al. built an 
all-sky network of absolute reference sources for spectro- 
photometry in the infrared. The tenth paper, C99, exposes at 
length their method: by fitting low-resolution composite spec- 
tra to infrared photometric measurements, they derive both ab- 
solute spectral irradiances and stellar angular diameters. Their 
grid stars were selected from the IRAS Point Source Catalog 
(hereafter IRAS PSC) by applying stringent quality criteria 
that qualify most of these stars as interferometric calibrators 
as well. Cat. 1 was then built by removing the very few stars 
that would depart from a UD at low resolution. 

In their original work, C99 used the following 6 criteria in 
order to select calibrator candidates in the IRAS PSC: 

1. The IRAS flux at 25 //m must be greater than 1 Jy; 

2. IRAS infrared colors are restricted to one particular quad- 
rant, corresponding to giant stars; 

3. IRAS candidates must be as normal as possible: sources 
flagged as non-stellar, variable, emission line or carbon 
stars are rejected; 

4. Total IRAS fluxes within a radius of 6' must be less than 
5 % of the IRAS candidate corresponding fluxes; 

5. IRAS stars must not be associated with a small extended 
source. The contribution of infrared cirrus must be less than 
5% of the candidates infrared flux; 

6. To reduce the chance of selecting variables, candidates 
must fall in one of the following spectral and luminosity 
classes: 

- A0-G9 and II-IV, 

- KO-MO and III-V. 

Walking in the tracks of C99, we applied the previous set of 
criteria to the IRAS PSC, modifying criterion 1 to F25 < 1 Jy 
in order to sort out new stars. We also required the stars to be 
identified in the SIMBAD database to make possible further 
quality filtering (criteria 3 and 6). A total of 8313 stars passed 
criteria 1 and 2, as well as the identification by SIMBAD. In the 
selection process, we found the 5th criterion - meant to remove 
stars blended with interstellar cirrus - severely impairs the sky 
coverage around the Galactic Equator Therefore, we decided 
to leave aside that criterion, and to rely instead on the quality 
of the spectrum fit to the infrared photometric measurements. 

4.3. Discarding candidates that depart from a UD 

Multiplicity and variability are the main causes of departure 
from a UD that one would like to avoid. Ideally, multiple stars 



should be rejected according to the separation and contrast be- 
tween the components. As thresholds on these quantities could 
only be determined according to the instrument configuration 
and characteristics (baseline, wavelength, accuracy, etc.), we 
chose to remove all known spectroscopic and astrometric bina- 
ries on the basis of their SIMBAD object type. However, be- 
cause of SIMBAD's hierarchical classification, a star can have 
a companion without being classified as a binary. For this rea- 
son, we also used the Catalog of Visual and Double Stars in 
Hipparcos (Dommanget and Nys 2000) to remove all stars with 
a companion within 2'. This distance represents a safe upper 
limit for any pointing devices or tip-tilt servos. In addition, 
the Hipparcos Catalog (Perrymanet al. 1997) was used in or- 
der to remove known variables, as well as stars with a positive 
"Proxy" flag, ie having field stars within 10". The 2' and 10" 
thresholds make the catalog useful for a wide variety of instru- 
ments without impacting significantly on the number of stars 
left in the catalog (only 12 % of the stars were discarded). 

In order to pin down more variable stars or unresolved bi- 
naries, we used several radial velocity catalogs, starting with 
the catalog by Malarodaet al. (2001) who collected all radial 
velocity measurements from the literature until 1998. Every 
entry of this catalog contains the nature of the object ac- 
cording to every author. If a single author suspected a star 
to be a spectroscopic binary, we discarded that candidate. 
Furthermore, we looked at de Medeiros and Mayor (1999) and 
Nidever et al. (2002) for post- 1998 measurements and used 
them in same manner as Malarodaet al. (2001). This process 
led to remove an additional 1 % of our candidates. 

A last cause of departure from the symmetric UD model 
would be oblateness due to fast rotation. This should not be a 
concern in our stellar sample as it contains only G8-M0 giants 
expected to be slow rotators. We checked that by using vsin(/) 
as a proxy for the rotation rate: indeed, for all 36 stars with 
measurements, we found that vsin(;) < lOkm/s. 

As a general remark, let us point out that considering the 
large number of candidate stars, it was out of question to check 
every one of them in the literature, and to perform dedicated ob- 
servations where it would have appeared necessary. Therefore, 
we would like to draw the attention of the reader to the fact that 
a few stars which are improper for calibration may remain in 
our published list. 

4.4. Photometry 

Deriving angular diameters with C99's method requires abso- 
lute photometric data between 1.2 and 35 yum. We used 12 fim 
and 25 fim photometric data from the IRAS PSC (our input 
catalog), as well as those from the IRAS FSC (Faint Source 
Catalog). We extracted the near-infrared photometry, H and Ks 
magnitudes, from the Two Micron All Sky Survey Point Source 
Catalog (hereafter 2MASS PSC) with the absolute cafibration 
by Cohen et al. (2003). The 2MASS PSC was the only near- 
infrared catalog that matches our needs in terms of sky cover- 
age. However, we had to deal with a poor photometric precision 
{ctk = 0.25) for the bright stars (K = 3^) we are interested in, 
as they tended to saturate 2MASS detectors. FinaUy, the MSX 
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infrared Astrometric Catalog (Egan and Price (1996)) provided 
6 photometric data points in the 4.22-4.36/im (called Bl band), 
4.24 - 4.45 //m (called B2 band), 6.0 - 10.9 jim (called A band), 
1 1 . 1 - 1 3 .2 yL/m (called C band), 1 3 .5 - 1 6 .9 yum (called D band) 
and 18.1 - 26.0 ;um (called E band) infrared bands. These last 
measurements were converted into absolute photometric mea- 
surements using Cohen et al. (2000) calibration. 

We looked for interstellar extinction by computing the 
spectral class photometry (especially B-V), comparing it to its 
expected value. If is the extinction in magnitude, the in- 
terstellar medium is characterized by Ry - AvI{Ab - Ay) - 
Av/Eb-v- We used an average of Ry - 3.1. Then, thanks to 
the tables of A^/Ay from Cardelli et al. (1989), we corrected 
the H and K magnitudes for reddening. However, this correc- 
tion always stays within the photometric error bars, because our 
candidates are bright stars that are close and for which 2MASS 
photometry is not very accurate. Therefore our relatively crude 
estimation of the reddening is appropriate. 

4.5. Deriving limb-darkened diameters from stellar 
templates 

Stellar templates computed by C99 consist of absolute compos- 
ite spectra for different spectral types and luminosity classes, 
along with limb-darkened (LD) angular diameters 0. Because 
C99 used plane-parallel Kurucz model atmospheres, their LD 
angular diameters {(p) correspond to Rosseland angular diame- 
ters, associated with a Rosseland optical depth of unity. Indeed, 
this diameter corresponds to the zone where photons from the 
continuum are emitted, and where the limb brightness drops to 
zero. In the case of plane-parallel models the Rosseland diam- 
eter is equal to the limb darkened diameter (see discussion in 
Wittkowski et al. (2004)). 

After retrieving the spectral types and luminosity classes 
of our candidate calibrators from SIMBAD, we fitted the 
absolute template corresponding to a specific star to its in- 
frared photometric measurements. All photometric measure- 
ments were converted into isophotal quantities, following the 
calibration prescribed by C99 for IRAS PSC and FSC, by 
Cohen et al. (2003) for 2MASS PSC and Cohen et al. (2000) 
forMSX. Fig. 4 shows two examples of stellar template fits. 

The scaling factor resulting from the fit, I, yielded the LD 
diameter of the star 

(p^ yflX 0,-ef . (6) 

For every diameter, we computed an associated error by taking 
into account the contributions of the fitting formal error, cr/, 
corresponding to - X^i^ ± cr/)l - 1' ^"d of the original 
template diameter error, cr^^^f, according to 




We could not find in C99 the angular diameter errors of 
their original templates. By comparing the diameters published 
by these authors and the ones computed by our implementa- 
tion of their method, this internal error was determined to be 
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Fig. 4. Examples of stellar template fits to photometric measurements 
converted into monochromatic irradiances at isophotal wavelengths. 
The data points stand for (from left to right): H, K, (2MASS), Bl, 
B2, A (MSX), Fn^m (IRAS PSC and FSC), C, D, (MSX) and Fzs^m 
(IRAS PSC and FSC). Left: HD 15248, a KlIII star Of? = 0.66). Right: 
HD 164724, a MOIII star (x'' = 0.83). 

(^(prdl'pKf - 1.03%. This value appeared to be independent of 
spectral type and luminosity class. 

As mentioned at the end of Sect. 4.2, we used the quality of 
the template fit as a filtering tool to ensure the consistency of 
our method. Namely, if the reduced;^^, i.e.xi - X^/iN^P+ 1)' 
where N is the number of data points and P the number of 
parameters, is greater than 1, we exclude the star from the final 
selection. 

4.6. From limb-darkened to uniform disk diameters 

We converted LD diameters into UD diameters in J, H and K 
bands by using the following four-coefficient limb-darkening 
law from Claret (2000) 

4 

/(y^)= 1- J]fl,(l-/^2), (8) 

A-=l 

where / is the specific intensity and fi the cosine of the angle 
between the line of sight and the perpendicular to the stellar 
surface. Claret (2000) has tabulated the coefficients a^, accord- 
ing to the stellar effective temperature T^ff and surface gravity 
log(g). 

Only 65 of our candidates could be found in the 
Cayrel de Strobel et al. (2001) catalog of effective temperature 
and surface gravity measurements, all of them with a spec- 
tral type G8-K4. For all others, we estimated T^ff by using the 
Bessell et al. (1998) polynomial relation based on the color in- 
dex V - K. When no measurements were available, we used 
the typical surface gravity for the spectral class, as well as 
solar abundances. This is legitimate since the infrared limb- 
darkening is not very sensitive to these parameters. 

Taking into account a center-to-limb variation, the visibility 
function becomes 

yhDW -J . (y) 

UD diameters in J, H and K, were obtained by fitting Vud 
(Eq. 3) to Vld (Eq. 9) in the domain of validity of the UD 
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Fig. 5. Histograms of magnitudes 



model, i.e. > 0.4. In this domain, the relative difference 
between the two visibility profiles is less than 10"'' (see Fig. 1 
for the most darkened disk of our sample, namely MOIII stars 
in J band). 



Fig. 6. Sky coverage of Cat. 2 in equatorial Hammer-Aitoff represen- 
tation. 



Dist. to the closest calib. 


Northern Hem. 


South. Hem. 


less than 3.75° 


45 % (15 %) 


85 % (53 %) 


less than 5° 


75 % (45 %) 


95 % (65 %) 


less than 10° 


99 % (85 %) 


100 % (89 %) 



5. Major characteristics of the catalog 

5.1. Electronic version 

An electronic version of the catalogue is available online from 
the Centre de Donnees astronomiques de Strasbourg (CDS). 
For every calibrator star, it provides: 

1. Identifiers: HD and HIP numbers; 

2. Coordinates: ICRS 2000.0 right ascension and declination, 
proper motion, and parallax; 

3. Physical properties: spectral type; 

4. Photometric measurements: B and V from SIMBAD, J, H 
and Ks from 2MASS PSC, 12 |/m from IRAS; 

5. Angular diameters: LD diameter and UD diameters for J, H 
and K. A single error is given for all of them, since the cor- 
rection on this quantity due to the center to limb variation 
is very small; 

6. Ancillary information: variability and multiplicity flags, 
variabiUty type from Hipparcos. This category is meant to 
distinguish between stars for which no variability nor com- 
panion has been detected and stars for which no informa- 
tion is available. 

5.2. Infrared magnitudes 

On Fig. 5, we plotted three histograms for J, H and Ks infrared 
magnitudes. The median magnitudes are 4.3 in J, 3.6 in H and 
3.5 in Ks. 

5.3. Sky coverage 

Cat. 2 overall sky coverage is such that there is always a 
calibrator closer than 10° whatever the location on the sky 
(Table 1). Moreover, if we keep only stars that can achieve a 
2 % visibility error for baseline greater than 200 m in K (145 m 
in H, 110 m in J), the calibration appears to be still feasible 
for more than half the sky. The northern hemisphere happens 
to be less populated than the southern one (Fig. 6), because 
more stars were filtered out in this part of the sky due to 



Table 1. Sky coverage characteristics: area of the sky having corre- 
sponding distance to the closest calibrator. Details for northern (posi- 
tive declination) and southern (negative declination) hemispheres are 
given. The number in brackets corresponds to a reduced version of 
Cat. 2, containing only stars that can achieve a 2 % visibility error for 
baselines greater than 200 m in K (which corresponds to 145 m in H, 
110 m in J). 



Band 


1 % 


2 % 




3 % 


J 


40 ^ 91m 


55 125 m 


65 


^ 150 m 


H 


55 ^ 121 m 


75 167 m 


85 


^ 200 m 


K 


75 ^ 163 m 


100 ^ 224 m 


120 


— > 266 m 



Table 2. Improvement of the median maximum baseline for J, H and 
K bands at the 1, 2 and 3 % visibility precision levels between Cat. 1 
and 2. 



an unsatisfactory spectral type identification with SIMBAD. 
For the southern hemisphere, spectral identification comes 
from the Michigan Catalogue of Two dimensional spectral 
types (Houk and Cowley 1975), which offers better spectral 
type identification than what is available from any studies in 
the northern hemisphere. 

5.4. Baseline range 

Fig. 7 shows a comparison between Cat. 1 and 2 in the calibra- 
tor decision diagram (observations in the K band with a 2 % 
precision on the visibility), and illustrates the significant gain 
due to this work. Cat. 2 can be used to achieve a 2 % visibility 
eiTor on the whole sky with baselines up to ~ 100 m, ~ 130 m, 
and ~ 180 m, in the J, H and K bands, respectively (see Table 2 
for details). Therefore, our extended calibrator catalog fulfills 
the needs for calibration of the VLTI, the CHARA array, and 
any other 200-m class interferometers in the near infrared. 
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10.0 



40m 




X = 2.2 |im ,7 

aWV2 = 2% 



are too bright for the major infrared surveys. From the decision 
diagram, it can be seen however that the gain to be expected is 
in a better sky coverage at 200-300 m baselines, rather than an 
extension of the maximum basehne for which Cat. 2 would be 
useful. 
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ae/e (%) 

Fig. 7. Comparison of Cat. 1 (gray dots) and 2 (black dots) in the cal- 
ibrator decision diagram. The 50, 100, 200 and 300 m curves plotted 
here, for the K band (A = 2.2/jm), corresponds respectively to 35, 75, 
145 and 220m in the H band (A = 1.6yum) and 25, 55, 1 10 and 165 m 
in the J band (A = 1.2 pm). 

6. Conclusion 

We have presented a list of 1320 stars selected to serve 
as calibrator stars for long baseline stellar interferometry. 
Among them, 374 come from the previous catalog compiled by 
Borde et al. (2002). The selection of the new stars was made 
according to the spectro-photometric criteria defined by C99, 
as well as additional ones that are specific to near-infrared in- 
terferometry. We used the work by C99 to estimate the angu- 
lar diameter from infrared photometric measurements extracted 
from IRAS, 2MASS and MSX Astrometric catalogs. Our new 
catalog of interferometric calibrators is an homogeneous set of 
stars, provided with limb-darkened angular diameters and uni- 
form disk angular diameters in J, H and K bands. Its intrin- 
sic characteristics - a small list of stars with accurate angular 
diameters and uniform sky coverage - inherited from the pre- 
vious version, were carefully maintained, while we included 
smaller stars in order to fulfill the needs of 200-m class inter- 
ferometers, like the VLTI in the H band (and redder) in the 
bright star regime (K :^ 8), and the CHARA array in the K 
band (and redder). In this paper, we improved the performance 
of Cat. 2 with respect to Cat. 1 by extending the selection to 
fainter stars (this corresponds to going downwards on the cali- 
brator decision diagram of Fig. 3). As all eligible stars from the 
IRAS survey have been used in Cat. 2, further extension is not 
possible in this direction. Improvements could still be obtained 
if/when more accurate photometry becomes available, which 
will result in smaller diameter errors (this corresponds to going 
leftwards on the calibrator decision diagram). For this, a ded- 
icated observing program is needed as most of these sources 
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